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Abstract 

We explore the dependence of the pT correlations in the event-by-event analysis of 
relativistic heavy-ion collisions at RHIC made recently by the PHENIX and STAR 
Collaborations. We point out that the observed scaling of strength of dynamical 
fluctuations with the inverse number of particles can be naturally explained by the 
formation of clusters. We argue that the large magnitude of the measured covariance 
implies that the clusters contain at least several particles. We also discuss whether 
the clusters may originate from jets. In addition, we provide numerical estimates of 
correlations coming from resonance decays and thermal clusters. 

Keywords: relativistic heavy-ion collisions, event-by-event fluctuations, particle 
correlations. 
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Recently new data on the event-by-event fluctuations have been provided by 
the PHENIX [1] and STAR [2,3] Collaborations, shedding more light on the 
previously accumulated knowledge in the field [4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19]. 
One of the most fascinating but intricate questions is whether the pt fluctua- 
tions in large windows of pseudorapidity and azimuthal angle at intermediate 
momenta can result from jets [20,21]. In this letter we explore the basic facts 
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of the recent data [1,3]. In particular, wc argue that since i) the mean and 
the variance of the inclusive momentum distribution are practically constant 
at low centrality parameters, then ii) the variance of the average momenta 
for the mixed events is practically equal to the variance of the inclusive dis- 
tribution divided by the average multiplicity. Moreover, and this is our basic 
observation, i) also results in the fact that iii) the difference of the experimen- 
tal and mixed-event variances of average pt, denoted as cr^y^, scales as inverse 
multiplicity, as seen in experiments [1,3]. A possible explanation of this scal- 
ing can be provided by clustering in the expansion velocity: matter expands in 
"lumped clusters" of chunks of matter, having close collective velocity within 
a cluster, which induces correlations. Moreover, we show that the value of a^yn 
is large at the expected scale provided by the variance of pt, which indicates 
that the clusters should contain at least several particles in order to combi- 
natorically enhance the magnitude to the observed level. We discuss whether 
jets may be responsible for the formation of the clusters. Finally, we compute 
numerically the value of (j^y^ coming from the resonance decays and from 
thermal clusters in statistical models of heavy-ion collisions. The found values 
of the covariance per pair are small, suggesting larger numbers of particles in 
clusters. 

We begin by exploring the PHENIX measurement [1] of the event-by-event 
fluctuations of the transverse momentum at ^/sNN = 130 GeV. To simplify 
our notation, the letter p is used to denote \pt\, Pi is the value of p for the ith 
particle, and M = YL^=iPi/n is the average transverse momentum in an event 
of multiplicity n. The PHENIX results are recalled in Table 1. Several features 
of the data are striking: the quantities (M) and Up are practically constant in 
the reported centrality range c = — 30%, 



(M) = const. , Gp — const, (at low c). (1) 

We call the range of c where (1) holds the "fiducial centrality range" - this 
is where our conclusions will be drawn. We note that for peripheral events 
incomplete thermalization can result in a different strength of px fluctuations 
[19]. Next, we observe that ctm — CFp/ \J {n) . More precisely, for the mixed 
events one finds the formula 



'M 



+ 



(2) 



working at the level of 1 — 2%. Finally, the difference of the experimental and 
mixed-event variances of average pt, denoted as cr^ym scales to a remarkable 
accuracy as the inverse multiplicity. 
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Now we proceed to elementary statistical considerations. Consider events of 
multiplicity (of charged particles) n and transverse momenta Pi,P2, ■ ■ ■ ,Pn- 
The multiplicity n and the momenta are varying randomly from event to 
event. The probability density of occurrence of a given momentum configu- 
ration is P{n)pn{pi, ■ ■ ■ ,Pn), where P(n) is the multiplicity distribution and 
Pn(pi, ■ ■ ■ , Pn) is the conditional probability distribution of occurrence of pi, . . . ,p,, 
in accepted events, provided we have the multiplicity n. Note that in general 
p depends functionally on n, which is indicated by the subscript. The normal- 
ization is 

^P(n) = l, dpi...dpnPn{Pl,---,Pn) = 'i-- (4) 

n •' 



The marginal probability densities are defined as 



Pn • • • ^Pn-k) = J dpn-k+1 ■ ■ ■ dpnpn{Pl, ■ ■ ■ ,Pn), (5) 

with k — 1, . . . ,n — 1. These are also normalized to unity, as follows from 
Eq. (4). Since the number of arguments distinguishes the marginal distribu- 

Table 1 

Analysis of the cvent-by-cvcnt fluctuations in the transverse momentum. Upper 
rows: the PHENIX experimental data at ^/sNN = 130 GeV [1]; middle rows: the 
mixed-event results; bottom rows: our way of looking at the data. One observes that 
to a good approximation cr^/™'^ — (Jp/{n) and u^y^ = (^if — Cm™'^) ~ Except 
for the flrst two rows, all values are given in McV. The errors in the last row reflect 
the unknown round-ofl' errors in the data of the upper and middle parts. 



centrality 


0-5% 


0-10% 


10-20% 


20-30% 


(n) 


59.6 


53.9 


36.6 


25.0 




10.8 


12.2 


10.2 


7.8 


(M) 


523 


523 


523 


520 




290 


290 


290 


289 




38.6 


41.1 


49.8 


61.1 




523 


523 


523 


520 


mix 


37.8 


40.3 


48.8 


60.0 


^P\j M + 


38.2 


40.5 


49.8 


60.8 




60.3 ± 1.6 


59.2 ± 1.5 


59.8 ± 1.2 


57.7 ± 1.1 
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tions ^\ in the following we drop the superscript {n — k). Further, we 
introduce the following definitions 

{p)n = J dppn{p)p, var„(p) = J dppn{p) {p- {p)nf , 

C0V„(pi,p2) = j dpidp2 (Pl - (p)„) (P2 - (p)n) Pn{Pl,P2)- (6) 

The subscript n indicates that the averaging is taken in samples of multiplicity 
n. We note in passing that the commonly used inclusive distributions are 
related to the marginal probability distributions in the following way: 



Pm[X) 



:^P{n) / dpi...dpn^S{x -pi)pnipi,...,Pn) ^^nP{n)pn{x), 

n •' i=l n 

/n 
dpi...dpn -Pi)^iy -Pj)Pn{Pl,- ■ ■ ,Pn) 

(7) 



which are normalized to (n) and (n(n — 1)), respectively. 
For the variable M — Yn^iPi/n we find immediately 



(M)=^P(n) / dp,...dpnMp^{p,,...,Pn)=T.Pi^)iP)n, 
n n 

{M^) = '£ P{n) jdpi... dp„M^p„{pi, . . . , p„) (8) 



n 



J2 C0Vn(Pi,Pj) + n(n - 1) (p) 



Next, we use the experimental fact that the variance of the momentum distri- 
bution and its mean are independent of centrality in the fiducial range, which 
allows us to replace the quantities (p)„ by (M) and {p^)n ~ {p)n t'Y the 
average multiplicity, denoted as . In this way we get 



2 2 P(n) ^ P(n 



J2 COV n{pi,Pj) 



(9) 



In the mixed events, by construction, particles are not correlated, hence the 
covariance term in Eq. (9) vanishes and 
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where in the last equahty we have used the fact that the distribution P{n) is 
narrow and expanded 1/n = l/[{n) + {n — (n))] to second order in [n — (n)). 
Comparison made in Table 1 shows that formula (10) works at the 1-2% 
level. In addition, since (Jp^i^n) is not altered by the event mixing procedure, 
subtracting (10) from (9) yields 



P(n) " 1 



^dyn = E^^ E COV„(pi,Pj) ~ -— COV{nM,Pj)- (H) 



Now we come to the physics discussion. The scaling (3) imposes severe con- 
straints on the physical nature of the covariance term. For instance, if all 
particles were correlated to each other, J27,j=i,j^i^'^'^ {n)iPi, Pj) would be pro- 
portional to the number of pairs, and a^yn would not depend on {n) at large 
multiplicities. A natural explanation of the scaling (3) comes from the cluster 
model, depicted in Fig. 1. The system is assumed to have clusters, each 
containing (on the average) r particles. Below we keep r = const, for simplic- 
ity. The particles are correlated if and only if they belong to the same cluster, 
where the covariance per pair is 2 cov*. The number of correlated pairs within 
a cluster is r(r — l)/2. Some particles may be unclustered, hence the ratio of 
clustered to all particles is {Nci)r/{n) = a. If all particles are clustered then 
a = 1. With these assumptions Eq. (11) becomes 




Fig. 1. The cluster model of correlations. Particles are grouped in A'^ci clusters, 
containing on the average r particles. The particles within a cluster move at very 
similar collective velocities, indicated by arrows. 
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which comphes to the scahng (3). An immediate conclusion here is that the 
ratio a cannot depend on (n) (in the fiducial centrality range) in order for the 
scaling to hold. 

The question now is whether we can use the above results to draw conclusions 
on effects of jets (minijets), which have been proposed as a possible explana- 
tion of the experimental data even at the considered soft momenta [20]. Jets, 
when fragmenting, lead to clusters in the momentum space. The resulting full 
covariance from jets is then A^ci,jetj(j — l)covY(n)^, where A^cijet is the num- 
ber of clusters originating from jets, j is the average number of particles in 
the cluster, and 2 cov-i is the average covariance per pair. The total number 
of particles produced from jets is Ai'djetJ- On the other hand, the commonly 
accepted estimate of the dependence of A^cijetJ on centrality is accounted for 
by the nuclear modification factor Raa multiplied by the number of binary 
nucleon-nucleon collisions A^bin- Since Raa depends on the ratio (n)/(n)pp, 
where (n)pp is the multiplicity in the proton-proton collisions, in a given pt 
bin one finds 



which complies to the scaling of Eq. (12). We stress that this scaling follows 
just from the presence of clusters, and is insensitive to the nature of their 
physical origin as long as one imposes N^x ~ (^)- In other words, as long as 
Eq. (13) is used, the explanation of the observed data in terms of quenched 
jets agrees with the cluster picture. However, the explanation of the centrality 
dependence of the pr fluctuations in terms of jets based solely on Eq. (13) is 
insufficient and not conclusive: any mechanism leading to clusters would do. 
Microscopic realistic estimates of the magnitude of cov-^ and j are necessary 
in that regard, including the interplay of jets and medium. For the current 
status of this program the user is referred to [21,22]. 

Before continuing the analysis of the cluster model in a more quantitative man- 
ner we need to consider the effects of acceptance and detector efficiency. This 
is particularly important in the event-by-event analysis, since the experiments 
select particles with very clearly identified tracks, and thus the detector effi- 
ciency, denoted by a, is small. The number of observed particles is proportional 
to a, and the number of pairs contributing to the covariance is proportional to 
a^. Thus Eq. (12) may be rewritten as o\,^ = ^ cov* = cov*, where 

"full" denotes all particles (that would be observed with 100% efficiency), 
while "obs" stands for the actually observed multiplicity of particles. Thus 




pp 



(13) 



cov* = a; 



2 (''^)obs 

'•^y°a(r- !)■ 



(14) 
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Our estimate for a in the PHENIX experiment is of the order of 10%, which 
together with the numbers of Table 1 gives 



, 0.035 GeV^ 

- • (15) 

In the considered problem the coefficient 0.035 GeV^ is not a small number 
when compared to the natural scale set by the variance a"^ ~ 0.08 GeV^. We 
recall that | cov* |< (7^. Comparing the numbers, we note that for r = 2 the 
value of cov* would assume almost a half of the maximum possible value. This 
is very unlikely, as argued in the dynamical estimates presented below, which 
give cov* at most 0.01 GeV^. Thus a natural explanation of the values in (15) 
is to take a significantly larger value of r. Of course, the higher value, the 
easier it is to satisfy (15) even with small values of cov*. We call this picture 
the "lumped clusters": lumps of matter move at some collective velocities, 
correlating the momenta of particles belonging to the same cluster, see Fig. 1. 

The above estimates were based on the PHENIX data [1] , however, very simi- 
lar quantitative conclusions can be reached from the recently published STAR 
data [3]. We note that the measure (ApiApj) used by STAR is just the esti- 
mator for cr^yn- Indeed, elementary steps lead to 



]\f 1 1 A^cvont /j2 

(Ap.Ap,> = - E (16) 

^ ' event ^ ' event f;—i ' fe 

Comparison to (9) leads immediately for a large number of events to (ApiApj) — 
(j^yjj. Now, taking the values of Table I of Ref. [3] and assuming a — 0.75 we 
find cov*(r - 1) = 0.058,0.043,0.035,0.014 GcV^ for = 200,130,62 

and 20 GcV, respectively. The value at 130 GcV is close to the value (15). 
Interestingly, we note a significant beam-energy dependence, with cov*(r — 1) 
increasing with ^^snn- This may be due to the increase of the covariance per 
correlated pair with the increasing energy, and/or an increase of the number 
of clustered particles. 

In the last part of this paper we present some dynamical estimates of cov* 
in thermal models. The first calculation concerns the role of resonances in px 
correlations. Clearly, a resonance, such as the p meson, decaying into daughter 
particles induces momentum correlations. We make a numerical calculation of 
this effect in the model of Ref. [23,24], using the formula 



(rp2 

^ !)■) 



COVr 



Id^Pl§^l'^Si^){p-p, 



(17) 
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where dNn/d^p is the resonance distribution in the momentum space (ob- 
tained from the Cooper- Frye formula as described in Ref. [25]), pi and p2 
are the momenta of the emitted particles, Ep is the energy of a particle with 
momentum p, and the function C represents the experimental cuts. We note 
that from now on the letter p, depending on the context, denotes the four- 
or three-momentum. The results of our numerical study show that for the 
resonance mass between 500 MeV and 1.2 GeV the covariance cov*^^ varies 
between 0.005 GeV^ at low masses to —0.015 GeV^ at high masses, chang- 
ing sign around 700-800 MeV, depending on the assumed experimental cuts. 
Thus, cancellations between contributions of various resonances are possible; 
in fact, a full-fledged simulation with Therminator [26] revealed a negligible 
contribution of resonances to the px correlations. Of course, the "lumpy" fea- 
ture of the expansion was not implemented in the calculation. Details of this 
study will be presented elsewhere. 

The second model of particle correlations assumes that the particle emission 
at the lowest scales occurs from local thermalized sources. Each element of 
the fluid moves with its collective velocity and emits particles with locally 
thermalized spectra. This picture was put forward as a mechanism creating 
correlations in the charge balance function [25,27] resulting from charge con- 
servation within the local source. Correlations between particles emitted from 
the same cluster come from the fact that those particles are emitted from a 
source with the same collective transverse velocity. The average number of 
particles r originating from such a local source determines the strength of the 
surviving dynamical fluctuation in the whole event, as discussed above. The 
covariance between particles i and j emitted from a cluster moving with a 
velocity u is 

_ j di:,u>^ j d^pM - ijp^))fr{p,)!d'pM - (p^))/r(P2) 



where /"(p) = (exp(p ■ u/T) ± 1)~^ is the boosted thermal distribution and 
dH^ denotes integration over the freeze-out hypersurface. The result turns out 
to depend strongly on the temperature. Considering the emission of corre- 
lated pion pairs one gets cov*(7r, tt) = 0.0034 GeV^ for freeze-out parameters 
corresponding to the single freeze-out model [23] (T = 165 MeV, average flow 
velocity 0.5c) and cov*(7r, tt) = 0.01 GeV^ for parameters corresponding to a 
late kinetic freeze-out (T = 100 MeV, average flow velocity 0.6c). For realistic 
values of thermal freeze-out parameters the experimentally estimated value of 
the covariance cannot be accounted for, unless the number of charged particles 
belonging to the same cluster is at least 4 — 10. 

In conclusion, we have found that in the fiducial centrality range the scaling 
of the cr^yjj for the pt correlations with inverse particle multiplicity indicates 
the cluster nature of the system formed in relativistic heavy-ion collisions. 
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The clusters may a priori originate from very different physics: jets, droplets 
of fluid formed in the explosive scenario of the collision, or other mechanisms 
leading to multiparticle correlations. A larger number of particles within a 
cluster helps to obtain the large measured value of cr? ^. 
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